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Abstract. We study a variational framework to compare shapes, modeled as Radon 
measures on M. N , in order to quantify how they differ from isometric copies. To 
this purpose we discuss some notions of weak deformations termed reformations 
as well as integral functionals having some kind of isometries as minimizers. The 
approach pursued is based on the notion of pointwise Lipschitz constant leading 
to a space metric framework. In particular, to compare general shapes, we study 
this reformation problem by using the notion of transport plan and of Wasserstein 
distances as in optimal mass transportation theory. 
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Introduction 

One of the main goal in shape analysis relies in detecting and quantifying differences 
between shapes. The interest for such studies concerns a wide range of applications, 
especially those within the computer vision community, in particular in pattern recog- 
nition, image segmentation, and face recognition (see [H]). In recent years many 
authors have focused their attention on the notions of shape space and shape metric 
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to the aim of establishing a general framework in which the analysis of shapes cru- 
cially depends on their invariance with respect to suitable geometric transformations 
(see [23l [57] ). A natural suggestion in this direction comes from continuum mechanics 
since the variational theory of elasticity can be used to compare the initial and final 
shape of a deformable material body, i.e. to establish how the two shapes differ from 
an isometry of the euclidean space. Therefore some authors begin to study the possible 
links between elastic energies and distances in shape spaces (see [23l [Mf [65]). 
On the other side, by arguing from a mechanical perspective, we know that a large 
class of physical manifestations (fractures, fragmentations, material instabilities) re- 
quire more general kinematical tools than those available in the context of Sobolev 
maps, hence it seems reasonable to exploit a more general mathematical framework to 
obtain more accurate descriptions of more complex physical problems. 
In this paper we model (material) shapes as Radon measures on subsets X,Y <Z M. N 
and study a variational model to the aim of quantifying how a target shape v on Y 
differs from an isometric copy of /i on X. To this purpose we scrutiny some notions of 
weak deformations, which we denote by the term reformations, as well as energy like 
(or cost) functionals having some kind of isometries between n and v as minimizers. 
In the first part of the paper (Sections 1,2,3) we study the variational problem of ref- 
ormation of two shapes fi and v through functions called reformation maps, while in 
the second part (Sections 4,5,6) we relax the problem by considering a formulation in 
terms of transport plans which leads to a variational framework as in optimal transport 
theory. 

For reader convenience we have added an appendix containing some basic tools from 
analysis in metric spaces. 

Description of the variational model and main results 

The general question addressed in this paper is to provide some variational tools 
able to quantify how two shapes X, Y C ~R N are close to be isometric. An usual 
way to compare the two shapes relies in considering Y = u(X) for maps u belonging 
to a suitable class of admissible maps. The two shapes are isometric if there exists 
u : X — > Y such that u(X) = Y and 

\u(x) — u(y)\ = \x — y\, Vx, y E X. 

Equivalently, the last condition means that the map u has bi-Lipschitz constant L — 1. 
Let us recall that a map u : x — > Y is said to be bi-Lipschitz with constant L if 

— \x — y\ < \u(x) — u{y)\ < L\x — y\, Vx, y E X. 

Therefore, the two shapes X, Y could be considered close to be isometric as the bi- 
Lipschitz constant L is close to one, so assuming the bi-Lipschitz constant as a quan- 
tifier of the closeness to the isometry. This approach has the disadvtange to involve a 
global condition. For instance, the shapes in Figure 5.2 looks very close to be isometric 
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but the bi-Lipschitz constant is quite large and far from L = 1, whatever the size of 
the bending part. To avoid this difficulty some localization procedure is needed. This 
can be done by an analytical approach. 

An isometry u is of course an affine map u(x) = Ax + b and Vw = A is an orthogonal 
matrix. Actually, under some regularity assumptions, by Liouville Rigidity Theorem 
the orthogonality of the Jacobian matrix characterizes the isometric maps (see also 
Theorem 13 .8p . Hence, a reasonable way to quantify how two shapes are isometric is 
that of measuring how Vw is close to be an orthogonal matrix. This program can be 
carried on by selecting a function W reaching its minimal value at the orthogonal ma- 
trices. Then, by Liouville Rigidity Theorem, it follows that the isometries characterize 
the minimal possible value of the functional I(u) = j n W(Vu) dx. This approach 
is pursued in [SS] for smooth 2-dimensional domains where the admissible maps are 
incompressible diffeomorphisms, i.e. u such that det (Vw) = 1. 

In order to characterize the isometries, a polyconvex function W having minimal 
value at orthogonal matrices is selected. Therefore, to quantify how two domains 
fli, f^2 C M. N are close to be isometric one considers the variational problem 



where T> denotes the set of incompressible diffeomorphisms. This approach has of 
course many restrictions. For instance, to compare a connected domain to a discon- 
nected one, or for non-smooth domains, many regularity questions arise. 
This discussion is reminiscent of the variational principles in continuum mechanics, 
indeed in the setting of nonlinear elasticity, in the undeformed state the body occupies 
a bounded open set Q C M. N , and one usually look for minimizers u : ft — > M. N of the 
stored energy I (it) = J n W(Vu) dx in an admissible class of deformations usually con- 
sisting of Sobolev functions which are locally orientation preserving, i.e. detVw(x) > 
for a.e. x G Q. 

A main goal of our approach relies in exploiting possible extensions of the variational 
scheme of elasticity (see also [28] for an analysis of this topic) in order to compare more 
general shapes as those in Figure 4.1, also allowing fragmentations. However, a purely 
measurable setting does not work to compare shapes as shown in Example 13.61 and, 
on the other hand, to compare a more extended class of shapes we have to reduce 
regularity requirements. So, a useful compromise relies in working on a general metric 
framework. 

We remark that an approach like the one followed in [SS] cannot be pursued in a met- 
ric framework, indeed the mapping A i— > (p(\\A\\) is polyconvex only if <p is a positive 
convex and strictly increasing function (see for instance [15J), therefore the minimal 
value cannot be reached at orthogonal matrices A, since they have ||y4|| = 1. 
We denote by V(X) the space of probability measures on X. Assume the material 
shapes are given by probability measures /i G V(X) and v G V(Y), (to fix ideas con- 




sider ii = £ N \_X, v = C N LY). 
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In this paper we assume the pointwise Lipschitz constant Lip («)(#) (see Definition II .ip 
as a local descriptor to measure how an admissible map u is an expansion or a contrac- 
tion. Note that the pointwise Lipschitz constant Lip(w)(x) coincides with the operator 
norm ||Vw(x)|| whenever u is a different iable map. Hence the local expansion and con- 
traction due to the map u at any point x are respectively represented by the functions 
e u {x) and c u {x) (see Definition 12. ip depending on Lip(w)(x). 
We require the admissible maps u : X —¥ Y satisfy the conditions 

u#^ = u } (0.1) 

e u <K, c u <H (0.2) 

for given H, K > 0. These maps will be called reformation maps, and the set of 
such functions will be denoted by Ref(/i; v) (see Definition 12. 2p . We consider the local 
reformation energy r u = e u + c u and the total reformation energy M{u) = f x r u (x) dfi, 
so in Theorem 13.121 we show that the variational problem 

minimize-f^w) | u G Ref(/i; u)} (0.3) 

admits solutions whenever Ref (/i; v) ^ 0. Therefore, to quantify how the two measures 
/i, v are isometric we look to the elastic reformation energy between /i and v defined 
by 

£(//, v) := mi{M(u) \ u G Ref(/i; u)}. (0.4) 

In Theorem 13.171 we show that the value of (10. 4p is attained if and only if the two 
shapes are isometric. 

In Section 4 we extend the scenario to deal with the case of non-existence of maps 
satisfying (12 .4p and this happens, for instance, when fragmentation occurs. In such 
a case the notion of transport plan coming from the optimal mass transportation or 
Monge-Kantorovich theory is useful. A transport plan between /i and v is a measure 
7 G V{X x Y), having /z, v as marginals, namely (7Ti)#7 = ji, (7r 2 )#7 = z/, where 7Ti j2 
are the projections of X x Y on its factors. The notion of transport plan could be 
considered as a generalization of a transport map, i.e. u : X — > Y, such that u#fi = z/, 
and so as a weak notion of reformation of [i into v. We shall refer to such measures as 
reformation plans. Actually, to every transport map corresponds the transport plan 
given by (/ x w)#/i. The shapes /z, v could be compared by looking the local mass 
transportation displacement in the target configuration. 

More precisely, by Disintegration Theorem (see [31 Section 2.5]) every transport plan 
7 G V(X x Y) can be written as 7 = f(x) <g) /i, where 

f:X^(V(Y),W) (0.5) 

is called disintegration map and V(Y) is equipped with the Wasserstein distance W. 
This point of view leads to formulate the reformation problem in terms of disintegra- 
tion map / and related metric expansion and contraction energies (see Definition I4.6P . 
So, in this setting, reformation maps take value in the space of probabilities, endowed 
with the Wasserstein metric, over the target domain. The main advantage of this 
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approach relies in its generality and in its connections with fertile topics as optimal 
mass transportation and geometric measure theory. However, many interesting open 
questions arise as the regularity needed on / to capture relevant geometrical and phys- 
ical properties of the shapes. 

In Section 5 we show several examples of shape reformations attainable by disintegra- 
tion maps but not attainable by any regular transport map. 

In Section 6 we study the main aspects of the variational problems of reformation in 
the enlarged context of generalized reformations, showing in Theorem 16.41 how isomet- 
ric measures can be characterized by means of the reformation energy. In Theorem 16.81 
we prove the existence of minimizing reformation plans for a constrained variational 
problem. 



1. The pointwise Lipschitz constant 

In this section we introduce the notion of pointwise Lipschitz constant and scrutiny 
some properties related to this tool since it will play a crucial role in this paper. Notice 
that The pointwise Lipschitz constant is a useful tool in theories concerning Sobolev 
spaces in a metric setting (see [121 HH1 ETJ)- 

Definition 1.1. Let (X,d x ), (Y,dy) be two metric spaces and let f : (X,dx) — >■ 
(Y,dy)- The pointwise Lipschitz constant Lip(/)(x ) of f at x G X is defined by 



y d Y (f(x),f(x )) . J -I J 'J. 

T . / -w \ i hmsup : — ; r it in is a non-isolated point ,_, lN 

Lip(/)(x ) := { x-^xo d x (x,x ) (1.1) 

if xq is an isolated point. 

It is readily seen that Lip(/) is a measurable function. 

Lemma 1.2. Let L > 0, X C W N a convex set and let f : X — >■ (Y, d) be a function 
such that Lip(/)(x) < L Vx G X . Then f is L-Lipschitz. 

Proof. Let k > L and a, b G X be fixed. We define the set 

S = {te [0,1] : d(f(a),f(a + t(b-a))) < kt\b - a\} . 

Since Lip(/)(a) < \i then S ^ 0. Let r = sup S. By continuity of /, which follows 
by the boundedness of Lip(/)(x)( see also [18]), we have that r G S. We claim that 
r — 1. Otherwise, since Lip(/)(o + r(b — a)) < k, we find i G (r, 1) such that 

d(f(a), f(a + t(b - a)) < d(f(a), f(a + r(b - a)) + d(f(a + r(b - a)), f(a + i(b - a)) < 

rk\b — a\ + (t — r)\b — a\k = ki\b — a|, 

so t G S, in contradiction with the maximality of r. Letting k — » L + we get the 
thesis. □ 
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A result similar to the previous lemma holds true for quasi-convex metric spaces X 
(see HBJ). 

The pointwise Lipschitz constant is also related to the notion of metric differential (see 
S2J SSI [50]). A function / : X C M N ->• (Y, d) is said to be metrically different iable 
at a point x G X if there exists a (unique) seminorm on WL N , denoted by MD(f,x ), 
such that for every y,z G X the following formula holds true 

d(f(y), f(z)) - MD(f, x )(y - z) = o{\\y - x Q \\ + \\z - x \\). (1.2) 

Let U C M. N be an open set and let / : U — > (Y, d) be a Lipschitz function. Hence, for 
every fixed p EY the function 

x ^ d(f(x),p) : U R+ (1.3) 

is a Lipschitz function and by Rademacher Theorem it is a.e. different iable in U. 
Moreover (see [121 113] )> it turns out that / is metrically differentiable at almost every 
point. 

The following lemma establishes a link between the pointwise Lipschitz constant, 
the distance function (see [6] for dual Banach spaces) and the metric differential (see 
(O)!. 

Lemma 1.3. Let f : U C M. N — > (Y, d) be a Lipschitz function over a separable metric 
space Y . Then, for a.e. xq G U it results 

Lip(/)(xo) = sup||W(/(x ),p)||. (1.4) 

Proof. In the above formula we assume that ||V(i(/(xo),p)|| = if the function x i— > 
d(f(x),p) is not differentiable at x . For p G Y we compute 

iTiMti \ \ \ v d (f( x o + tv),p) - d{f{x ),p) 

(Vd(f(x ),p),v) = hm < 

t->o+ t 

< hm d{f{Xo + tv)J{Xo)) <Lip(f)M. 
t^o+ t 

Taking the supremum with respect to v and then respect to p, the inequality 

Bup||Vd(/(z ),p)ll <Lip(/)(z ) 

follows. To get the opposite inequality, we use a slight modification of the proof of [TJ 
Theorem 4.1.6]. Since Y is separable, we fix a countable dense set {p n } C Y, then for 
every x\,X2&U we have 

d(f(xi), f(x 2 )) = sup |d(/(a?i),p n ) - d(f(x 2 ),p n )\- (1.5) 

n 

Consider the Lipschitz function <p n (t) = d(f(xo + tv),p n ) and set m{t) = sup n |0 n (t)|. 
Observe that \m(t)\ < Lip(/). By the Lipschitz condition, we may suppose that, for 
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every n G N, t = is a differentiability point for ip n and also t = is a Lebesgue point 
for m G By (11. 5p we obtain 

d(f(x + tv),f(x )) 1 /■* 1 /■* 

< sup - / |0 n (sJ|as < - / m(s)ds. 



t 

Letting t — )• + we get (see Prop. 1 and Th. 2 of [43J) 

MD(f,x )(v) < m(0) < sup ||Vd(/(x ),p„)||. (1.6) 

n 

On the other hand, by (ll.2p we get 

d(f(x), f(x )) _ MD / x-xp \ + o(\\x - x \\) 



\\X — Xq\\ \\\x — x \\ / \\x — X Q \ 

Letting x — > x , by (II. 6p we get 

Lip(/)(x ) < sup \\Vd(f(x ),p n )\\- 



□ 



Observe that by the proof of the previous lemma, the following equality also holds 
true 

Lip(/)0r o )= sup \MD(f,x )(v)\. (1.7) 

r, H=l 



Lemma 1.4. Let (Y,d) be a separable metric space. Assume U C ~R N is an open 
se t, (/n)neN be a sequence of (locally) equi- Lip schitz functions /„:[/—>■ (Y, d) and let 
f : U — > (y, d) . If f n —> f (locally) uniformly on U then 



Lip(/)(x) dx < liminf / Lip(/ n )(x) dx. 



u 

u 



Proof. By uniform convergence / is a (locally) Lipschitz function, moreover we have 
that d(f n (-),p) — ^ d(f(-),p) weakly* in W^iJJ). Therefore, for every p, by weak l.s.c. 
of the gradient norm (see also Ch. Ill Th. 3.3 of [59J) and using (ll.4p we get 

\\Vd(f(x),p)\\dx< liminf / \\Vd(f n {x),p)\\dx < liminf / Lip(/ n )(x) dx. 

n— >+oo Jjj ri— >+oo Jjj 

Since Y is separable, as in the proof of Lemma fl~3| denoting by g n = || Vd(f(x),p n )\\, 
by (ll.4p we may assume that Lip(/)(x) = \im n g n (x). Moreover, observe that 1^(^)1 < 
Lip(/). Hence, passing to the limit under the integral sign we finally obtain 

/ Lip(/)(x) dx = lim / g n (x)dx < liminf / Lip(/ n )(x) dx. 

Jjj n^ + OO Jjj H^ + OO Jjj 

□ 
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Remark 1.5. The above Lemma holds true of course for the function Lip p (/), for any 
p > 1. If Y C M. N , the uniform convergence can be replaced by the weak convergence 
on the Sobolev space W 1,P (U). In such a case, Lemma rOl just states the lower 
semicontinuity of the p-Dirichlet energy in Sobolev spaces, since if u : R — > ~&L N 
is differentiable at x then hvp(u)(x) = \\Vu(x)\\, see also [59j Ch. 3 Theorem 3.3]. For 
a related semicontinuity result see also [60]. Notice that if U C M. N is replaced by a 
general metric space (X, dx) the above arguments cannot be applied. In this case of 
generality we are forced to deal with metric Sobolev spaces (see Appendix IB]) . 

2. Reformation maps 

In this section we introduce the class of reformation maps and establish some prop- 
erties of these functions. Though the definition of reformation map holds for general 
metric measure spaces, as a first step we restrict our analysis to the euclidean frame- 
work of subsets of R . 

To fix ideas consider Q C M. N an open bounded connected set, X = Q, Y C M. N . 
Moreover, /i, v are given Radon measures on X and Y respectively, with normalized 
masses, so that fi £ V(X) and v £ V(Y). 

Definition 2.1 (Expansion and contraction energy). Let Xq £ X and u : X — > Y . 

The pointwise expansion energy of u at x is defined by 

\u(x) — u(xq)\ 

^ u \Xq) '■= Lip(w)(£ ) = hmsup : : . (2.1) 

x— >xq | X Xq I 

The pointwise contraction energy of u at x is defined by 

c u {xo) ■= limsup — j^r — X °, . . . (2.2) 

x^-xq \U{X) — U{Xq) I 

The pointwise reformation energy of u at Xq is defined by 

r u (x ) = e u (x ) + c u (x ). (2.3) 

Definition 2.2 (Reformation maps). We shall call reformation map any map u : X — > 
Y such that the following conditions hold true: 

u#n = is, (2.4) 

Vx £ X 3H,K,r > 0s.t.e u {y) < K, c u (y) < H £ B{x, r) n X. (2.5) 
We shall denote by Ref(/z; v) the set of reformation maps between fi and v. 

The point in the above definition is that the constants H, K may depend just on 
the point x an not by the map u £ Ref(/i;z/) or the radius r > 0. Therefore, the 
reformation maps are characterized by locally bounded expansion and contraction. 
Notice that, by the bounds ( 12. 5p . any u £ Ref(/x; v) is continuous and, by Stepanov 
Theorem (see [21]), is a.e. differentiable in Q. In particular, by Lemma [TT2l reformation 
maps are locally Lipschitz on Q. 
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Remark 2.3. In a mechanical perspective, the constraints stated in (12. 5ft could be 
considered as a bound on the maximum expansion or contraction experienced by the 
material fl. In this setting, the assumption that of the constants H, K do not depend 
on the map u in (12. 5p corresponds to a constitutive property of the material under 
consideration. We point out that some bounds as in (12.51) are in some sense necessary 
to control the geometry of the reformations. For instance, in the case of v = 5 yo 
we have e u = 0, c u = +oo for any map u satisfying (12. 4p . On the other hand, 
mapping a bar into a bended one (see Fig. 2.1) by two piecewise isometries u%, u 2 
such that Ui(x ) ^ ^2(^0), we necessarily create a fracture at the point x . It results 
e«(x ) = +00 at the discontinuity point x . See also Example 13.61 Therefore, roughly 
speaking, the bound c u < H means no collapsing, while e u < K means no fractures. 




Figure 2.1. Mapping a bar into a bended one. 

Remark 2.4. The constraint c u < H in (12. 5 h is related to inversion properties, both 
local or global, of reformation maps, see [2U |25j HQ]. Observe that for different iable 
maps with non- vanishing Jacobian we always have 

c u = ||(Vm) _1 ||, e u = \\Vu\\. 

Therefore, it would be an interesting question to consider pointwise constants H, K 
depending also on the map u. This is similar in spirit to the passage from functions 
with bounded distortion to functions with finite distortion (see the monograph [39J). 
However, in such a case, inversion properties becomes more subtle and further as- 
sumptions are needed, see for instance [MJ H7J [58] for inversion results of Sobolev 
maps. Anyway, the main interest of reformation maps relies in this perspective in 
considering just metric objects. See Section 4 for an extension in a metric setting. 
However, in a purely metric framework, such pointwise conditions are not enough to 
guarantee inversion properties. Consider for instance the map u : R — > M., u(x) = \x\ 
having e u = c u = 1 at every point. Therefore, also uniform bounds alone are not 
enough to get satisfactory inversion properties. Under differentiability assumptions in 
R , pointwise bounds are in fact enough, see Lemma [2.91 However, in general this is 
not true. For instance, see |54j, it is possible to find everywhere different iable maps 
with everywhere invertible differential on Hilbert spaces which are not neither open 
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or locally one-to-one. In the metric setting different restrictions arise, see [25] for a 
detailed discussion. 

Remark 2.5. The uniform bounds (12. 5 j) are also related to quasi-isometries, see for 
instance [HI ED] • In such case uniform bounds 

m < D-f(x) < D + f(x) < M, 

where 

n -ff \ r ■ r \f(x) -f(x )\ +f v \f(x)-f(x Q )\ 
D j(xo) = hmmf j j , D^j(x ) = hmsup j j 

x-^xo \X — X \ x ^ Xo \X — Xq\ 

denotes the local distortion of distances, are required. Observe that e u (x) = D + f(x), 
while c u (x) = D -)( XQ y Indeed, 

| a; — x \ I ^ 1 



\f(x) - f(x Q )\ i/w-/(»o)i - inf \m-H*o)\ ■ 

\x-x \ \x—xo\ 

Taking the supremum over B = B(xq,t) and then letting r — > + we get c u (x) < 
D -j( XQ ^ ■ In a similar way the opposite inequality follows. 

Remark 2.6. The mass conservation property (12.41) is a generalized version of incom- 
pressibility and it can be always satisfied (provided /i has no atom, see [53]) by some 
measurable map u. Actually, condition ( 12. 4ft is equivalent to the following change of 
variable formula 



f(u(x)) dfx = / f(y) du, (2.6) 
x Jy 

for every continuous function / : Y — > M.. 

Since Ref (fi; v) is made by nice functions, to compare the present approach with 

the classical ones (see also [65J ) , one should also require the condition u(X) = Y in 

the case of fi — £ N \_X, v = C N L.Y . However, we point out that this surjection 

requirement is actually a severe constraint. Indeed, for instance, to find Lipschitz 

functions u : X — Y Y, u(X) = Y for general compact sets in dimension N > 3 (see 

[Tj), as far as we know, is still an open question. Moreover, also bi-Lipschitz functions 

between nice sets are not easy to find (see [20l [31]). Anyway, in such a case many 

restrictions on the target space Y may be needed (connectedness, for instance). 

In the following we prove some properties enjoyed by reformation maps. A first 
estimate easy to verify is an immediate consequence of Definition 12.21 and is given by 
the following proposition. 

Proposition 2.7. Let u G Ref(/i; v). Then, for every xq G X there exists r > such 
that 

— \x — xo\ < \u(x) — u(xo)\ < K\x — xq\ VxgX D B(x , r). (2.7) 
H 
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Lemma 2.8. Let u G Ref(/i;z/). Then u is a discrete map, i.e. for every y G Y , 
u~ 1 (y) is a finite set. 

Proof. Let y G Y. By ( 12. 7p . if x G u~ 1 (y) we have that u(x) ^ u(x ) for every 
x G B(xq, r). Hence, Xo is an isolated point of u _1 (?/). We claim that is a finite 

set. Indeed, otherwise, since X is compact, we find a sequence x n — > x G X such that 
x n G By continuity of u we have also x G This is a contradiction 

since x is an isolated point of u~ 1 (y). □ 

For discrete continuous maps the local degree or local index i(xo,u) of u : X — > Y 
at xo G X is defined as follows 

i(x , u) = degO(xo), u, B(x , r)), (2.8) 

where deg(y, m, 5) denotes the topological degree (we refer to [HI [59] for an introduc- 
tion to degree theory). 

Let u^iy) = {x\, ■ ■ ■ ,Xh}, we have the following relation 

h 

deg(y,u,Y) = J2^j^). (2.9) 

3=1 

Observe that if u is locally injective in a neighborhood of x G X then |i(x,w)| = 1. 
We say that u is a sense-preserving (reversing) continuous map if the local index i(x, u) 
has constant sign in X. Notice that each homeomorphism on a domain is either sense- 
preserving or sense-reversing (see [191 Theorem 3.35]). Moreover, sense-preserving or 
sense-reversing different iable maps have constant Jacobian sign, (see [19j Lemma 5.9]), 
since 

i(xo,u) = sign(Ju(xo)), (2-10) 
where Ju := detVw, providing Ju(xq) ^ 0. 



Lemma 2.9. Let u G Ref(/i; v). If u is differentiable at x$ G Q then Ju(xq) ^ 0. 

Proof. Suppose by contradiction that Ju(xq) = 0. Then we find a vector \v\ = 1 such 
that Vw(xo) • v = 0. Fixed e > 0, since w is differentiable, there exists 5 > such that 
|u(xo + tv) — u(xq)\ < et, whenever \t\ < 5. On the other hand, by (12. 7p . there exists 
< t < 5 such that < \u(x + tv) — u(x ) \ < et. Hence < e. Letting e — > + we 
get a contradiction. □ 

By the previous lemma, any u G C 1 (X;F) fl Ref(/i;^) (or even an everywhere 
differentiable reformation map) is locally invertible on Q (see [32] for a related inversion 
result and [51] for an elementary analytical proof). 

If u is only a.e. differentiable, by Lemma [2.91 we have Ju(x) ^ for a.e. x. However, 
it is well known that in general this condition does not ensure the local invertibility 
of Sobolev maps (see for instance [46]). By the way, the condition Ju > on an 
open set fl is a standard requirement (see for instance [58]), ensuring that u is locally 



12 



LUCA GRANIERI, FRANCESCO MADDALENA 



invertible for a.e. x G VL. The restriction to sense-preserving maps is also made in [25] 
to derive local inversion problem. To this aim, also assumptions on the boundedness 
of HK < M, for sufficiently small M are necessary. In our context, since we are 
interested in comparing domains, also in a metric framework, restrictions to open 
maps seem more natural. 

For reader convenience we recall the following well known result 

Theorem 2.10. Let u : Q — > W N be a continuous open and discrete map. Then u is 
sense-preserving or sense-reversing. 

Proof. We essentially follow (8]. First, by Cernavskii Theorem (see also [32]) it follows 
that dim{u(B u )) < N — 2 and dim(B u ) < N — 2, where B u denotes the branch set 
of u, i.e. the set of points x G Q such that u is not a local homeomorphism at x, 
while the dimension dim is understood in the sense of [37J. Here, of course, we assume 
N > 2. For N — 1 the result is not true as it happens for instance for u(x) = \x\. 
Since dim(Q) = N, it follows (see [37J Theorem IV]) that Q \ B u is connected. Since 
u is a local homeomorphism on Q \ B u we have 

i(x,u) = c \/xett\B u , (2.11) 

where either c = 1 or c = - 1. Let xq G B u and yo = u(xq). By ( 12. 7p we find a 
closed ball C := B(xo,r) such that u(x) ^ u(xq) Vs G B(xo,r) \ {x }. Also, in a 
neighborhood V of ?/o we have 

deg(y, u, C) = i(x , u) My G V. 

Consider now U = u^ 1 ^) which is open in Q. Since u is open, it follows that u(U) 
is also open in u(Q). Since dim(u(B u )) < N — 2 we easily find y* G u(U) \ u(B u ). 
Let u~ l {y*) PI C = {xi, . . . , cfl \ Therefore, by (12. lip m) = c, for every 
i — 1, . . . , k and, since y* G V, by (I2.9p we compute 

k 

i(x ,u) = deg(y*,u,C) = S ^i{x i ,u) = kc. 

i=l 

The above formula shows that the index i(x,u) has constant sign. □ 

Remark 2.11. The previous theorem holds true in the case X, Y are topological man- 
ifolds. 

The following statement employs the same arguments of [59, Ch. II Theorem 6.6 ]. 

Lemma 2.12. Let u : Q — > M. be a discrete sense-preserving (reversing) continuous 
map such that \i(xo,u) \ = 1. Then u is infective in a neighborhood of Xq. 

Proof. Suppose i(xo,u) = 1. The other case is analogous. By contradiction, if u is not 
injective we have two distinct sequences (x^) ng N, (x^) ng N converging to x such that 
for every n G N: u(xl) = u(x^) = y n . By continuity we also have u{x^) — > yo = u(xq). 



ELASTIC REFORMATIONS 13 

Since the degree is constant in a neighborhood of yo, for n G N large enough and 
suitable small radius r > we have 

deg(y n ,u,B(x 0l r)) = deg(y , u, B(x , r)) = i(x 0l u) = 1. (2.12) 

On the other hand, since u is sense-preserving and by (12.91) we have 

deg(|/ n , u, B(x , r)) > i(a£, u) + i(x 2 n , u) = 2, 

contradicting (12.121) . □ 

We get the following inversion properties. 

Theorem 2.13 (Incompressible maps are invertible). Assume u G Ref(/i;^) is an 
incompressible map (see for instance [65],), i.e. 

\Ju(x)\ = 1 fora.e. x G O. (2.13) 

Then u\n is globally invertible. 

Proof. Of course, condition (12.41) holds true for injective such maps u. Anyway, let us 
begin by showing that u is an open map. Let U be an open subset of Q. We have to 
prove that V = u{U) is open. Let yo G V, yo = /(^o) for a xo G U. By (12. 7p we find 
C = B(xo,r) such that yo ^ u(dC). Therefore, it is well defined deg(yo,u,C) which 
is constant in a neighborhood of y . If u is different iable at x , by Lemma 12.81 and 
Lemma [2.91 we have 

|deg(?/,M,C)| = \deg(y ,u,C) \ = 1, 

in a neighborhood V yo of yo. Since deg(y,u,C) ^ 0, it follows that V yo C w(C) C V, 
since otherwise the degree would be zero. 

Observe that since Xb(x) < Xu(b)(u(x)), by Lusin Theorem and (I2.6p . we have 

C N (B)= [ X B(x)dx< [ X u(B)(u(x))dx= [ X u(B)u(y)dy = C N (u(B)). (2.14) 
Jx Jx Jy 

On the other hand, denoting by N(y,u,K) := card{« _1 (y) D K} the multiplicity 
function, by the Area Formula we compute 

C N (u(B(x ,r)))< [ N(y,u,B(xo,r))dy= [ \Ju\ dy = C N (B(x ,r)). 

Jr n JB(x ,r) 

Therefore, for a.e. y G u(B(xo, r)), it results N(y, u, B(xo, r)) = 1. Hence, if xo is 
a non-differentiability point, by the Lusin (N)-property, there exists a differentiability 
point x G B(x , 8) of u and an open neighborhood V yo of y such that m(s) G V^ and 
deg(u(x), u, C) 7^ 0. As above it follows that V yo C u(C) C V. Hence wiq is open. By 
Theorem 12.101 we have that u is sense-preserving, or sense-reversing. Since actually it 
results |z"(zo,w)| = 1, the statement follows by Lemma 12. 12[ □ 



14 LUCA GRANIERI, FRANCESCO MADDALENA 

Remark 2.14. For reformation maps satisfying ( 12.131) . it also results = 1 for 

a.e. x G Q. Therefore, by Lemma 12.121 such maps u are a.e. locally invertible on 
fl As previously discussed, actually to obtain this invertibility property it is enough 
to require \ Ju\ < 1 a.e. in Q. Since \Ju\ < ||Vw|| , this happens for instance for 
reformation maps u satisfying the condition e u < 1. 

Theorem 2.15 (Small reformations are invertible). Let u G Ref(/x;z/) be such that 
e u < \pl. Then u is globally invertible. 

Proof. By the constraint c u < H we find (see [511 Prop. 1.1, Sec. 3]) an open dense 
subset U C X on which u is locally bi-Lipschitz. It follows that the multiplicity 
function N(y,u,U) is locally constant (see [2]). We first prove that u is globally 
invertible on U. For y = u(x) G u(U), we have to prove that N(y,u, U) = 1. Observe 
that by the Domain Invariance Theorem, um is open. Let B = B(x,r) be a ball on 
which u is bi-Lipschitz. We may suppose that D = N(y,u,B) is constant on u(B). 
Since | Ju\ < ||Vu|| , by the Area formula we compute 

DC N (u(B))= [ N(y,u,B)dy= [ \Ju(x)\dx < [ \\Vu\\ N dx= [ e u {x) N dx. 
Jr n J b Jb J b 

(2.15) 

Observe that by the push-forward condition, as in (I2.14p . we obtain C N (B) < C N (u(B)). 
Therefore, if the map u satisfies the following small expansion condition 



e u (x) N dx < 2C N (B), (2.16) 

B 



then by fF2TT6D and (|2TT5]I we get 

DC N (u(B)) < 2C N (B) < 2C N (u(B)), 

hence the map u is globally invertible on U. By uniform continuity, u uniquely extends 
to the whole X and therefore letting global invertibility on A. □ 

Remark 2.16. By Area Formula, for B = B(x ,r) using (12.151) it results 

, \ m i £ N (u(B)) 

e x) dx > — — 

B u{) ~ C N {B) " 

By Lebesgue Differentiation Theorem this yields the a.e. inequality 

eu{x ) > ^hm ^ j := fi u {x Q ) » . 

The volume derivative fi u (xo) is also related to quasiconformal maps. We refer to 
BUSS, see also 



Theorem 2.17. [251 Th. II] Let u G Ref(/x;z/) be an open map such that HK < 2. 
Then u\q is locally invertible. 
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Proof. By Lemma 12.81 and Theorem 12.101 it follows that u is sense-preserving or sense- 
reversing. Hence, [23 Th. II] applies. □ 

The condition e u c u < M may be required to hold just a.e. by reducing the upper 
bound M < \/2 (see [25]). Observe that the results of [25] hold for dimension N > 2. 
Anyway, the map u(x) = \x\ in one dimension is not a counterexample to Theorem 
12.151 since u is not mass preserving around xo = 0. If Q is a ball, or in some classes 
of convex sets, for sufficiently small M the map u is actually globally invertible (see 
[261127]). 



3. The variational problem of elastic reformation 

Definition 3.1. Let u : X —*Y, [i G V(X), v G V(Y), such that u#fx = v . We define 
the total reformation energy M(u) of a reformation map u of \i into v as follows 

M{u) := / r u (x) dfi. (3.1) 
Jx 

To ensure that M(u) < +oo for every u G Ref(/i; u), we will always assume 

H(x), K(x)eL\X,fi), (3.2) 
where H, K are given in Definition 12.21 We have the following 
Lemma 3.2. 

c u {x) > —-^ Vx G X. (3.3) 

Proof. It suffices to recall that c u (x) = D - l f^ , see Remark [2751 □ 

Observe that for every u : X — > Y such that u#/j, = v, by Lemma 13.21 it results 
M{u) > 2. 

Actually, Definition 13.11 is motivated by the trivial fact that the real function f(x) = 
x + 1/x reaches its minimum value at /(l) = 2. Moreover, observing that at any 
x G X 

r u (x ) > e u (x ) + —— > 2, (3.4) 
e u {xo) 

we have that r M (x ) reaches its minimum value if u : X — > Y is an isometric mapping, 
i.e. \u(x) — u(y)\ — \x — y\, Vx,y G X. Therefore M{u) can be viewed as a measure 
detecting how u is far from being an isometric map. 
A useful property is the following 

Lemma 3.3. Let u : X — >■ V be a local homeomorphism. Then 

c u (x) = e u -i(u(x)) Vx G X. (3.5) 
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Proof. Fix x G X, 5i > and let B\ = B(u(x),5i). By the local homeomorphism 
condition, there exists a 5 > such that u(B$) C -Bi and it is invertible on 5^ = 
B(x, 5). For every y E B$ we have 

|y — a; I |M _1 ( , u(y)) — m _1 (-u(x))| Iw -1 ^) — -u _1 (-u(x))| 

< sup 



\u(y)-u(x)\ \u(y)-u(x)\ zeBl \z - u{x)\ 

Taking the supremum with respect to y G B$ and letting 5\ —> + , we get c u (x) < 
e u -i(u(x)). Analogously we deduce the opposite inequality. □ 

Definition 3.4. We define the elastic reformation energy between [i and v as 

£{ji, v) := mf{&(u) | u G Ref(/i; v)}. (3.6) 

Since the Monge transport problem is not symmetric, in general also the above 
elastic reformation energy is not symmetric. For instance, transporting a ball, say, 
into a Dirac delta we have v) = +oo. Reversing the shapes, we see that £ (v, //) 
has no meaning simply because Ref (v; fx) = 0. Moreover, also in nice cases, the matter 
is that transport maps could be not invertible. Assuming invertibility for u G Ref (//; u), 
setting v := u^ 1 , by using Lemma [3.31 we have 

M{u) — I e u d[i+ / c u dfi = / e u {u~ l {u{x))) dfx + / e u -i(u(x)) dfi — 
J x Jx Jx Jx 



e v -i(v(y)) dv + / e v (y) dv = / c v (y) dv + / e v (y) dv = M{v). 
y Jy Jy Jy 

Since v G Ref(i/;/i), we get £{\i,v) = £(v,fi). Therefore, symmetry issues essentially 

correspond to invertibility of maps. 

The question is now to establish conditions in order the infimum in (I3.6P is attained. 

We also want to prevent pathological situations as the one described in Example 13.61 

below in which the map u : X — > Y is merely a.e. continuous (it is actually a.e. 

invertible and differentiable). 

In order to select isometric maps through a variational property, we need to charac- 
terize the maps u : X — )■ Y such that r u (xo) = 2. 

Lemma 3.5. Let xq G X , u : X — > Y . Then r u (x ) = 2 if and only if 

We > : - — \x — xo\ < \u(x) — u(xo)\ < (1 + e)\x — xo|, Vx G X fl B(xo, r £ ). 
Proof. Assume r u (xo) = 2, then 

2 = e u (x ) + c u (x ) > e u (x ) H ^— > 2, 

so 

e u (xo) H 7 — r = 2 =>• (e u (x ) - l) 2 = =^ e u (x ) = c u (x ) = 1. 
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Figure 3.1. A piece-wise isometric map for the circle into a square. 



Fix e > 0, then e u (x ) < 1 + e implies that u satisfies 

\u{x) -u(x )\ < (l + e)\x-x \, Vx G XnB(x ,r £ ). (3.7) 

By using the condition c u (x ) < 1 + e, eventually by decreasing the radius r £ , we get 
the opposite inequality. Vice versa, if both the inequalities locally hold, then it results 
2 < r u (x ) = e u (x ) + c u (x ) < 1 + 1 = 2. □ 

Therefore, the maps u : X — > Y such that r u = 2 are in some sense pointwise locally 
quasi-isometric, (see [52] for the relation with quasi-conformal maps). 
In the following we shall try to characterize in a more precise way the reformation 
maps u G Ref(/i; u), if any, realizing the minimum energy level &(u) = 2. It is easily 
seen that 

£%{u) = 2 if and only if r u (x) = 2 for — a.e. x e X. (3.8) 

In particular, by ( 13. 4p (see also the proof of Lemma 133]) it results e u (x) = 1 for /i-a.e. 
x G X. Moreover, e u (x) < +oo implies u continuous at x. Then these reformation 
maps u are at least a.e. continuous functions. However, this mild regularity is too 
poor to preserve geometric (or physical) properties as we show in the next example. 

Example 3.6. Let Q C M. N be a smooth bounded open set and Q C R N be a cube such 
that C N (Q) = C N (Q), see figure 3.1. For n > 1 large enough, Q contains a certain 
number of disjoint squares Q n of length -. Then consider the map u n which moves 
by an isometry every square Q n inside Q in a disjoint way. On the remainder of Q, 
consider the contained squares Q m , m > n, and then the map u m which coincides 
with u n on the squares Q n and moves by an isometry the squares Q m inside Q in a 
disjoint way. By this procedure it is then defined a sequence (u n ) ne N- Taking the limit 
u = lim n ^ +00 u n we obtain a measurable map u : Q — > Q such that = u, where 
ji = C N \_Q, v = C N \_Q, and r u (x ) = 2 for a.e. x G Q. Therefore, every bounded 
smooth open set can be reformed into a square at minimal energy. 
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In order to preserve geometric and physical properties of the shapes under consid- 
eration, we then need more regularity on the admissible reformations. We have the 
following 

Lemma 3.7. Let xq G Q, u : X — > Y . If u is differentiable at Xq then 

r u (x ) = 2 Vu(x ) G O(N). 

Proof. By (13 .4p we have c u (x ) = e u (x ) = 1. Hence, for every v G WL N , taking 
x = Xq + 5v we get 

c u (x ) = e u (x ) = 1 =► \ Vu ^- v \ = i Vn(.„) G O(N). 

\v\ 

□ 

By Liouville Theorem (see for instance [13]) it follows that every u G C 1 (X; Y) such 
that M(u) = 2 is actually an isometry. There are several generalizations of Liouville 
Rigidity Theorem, however (see [131 HTJ [22]) these results are not directly applicable 
in our context since they generally require a constant sign for the Jacobian, as the 
condition Vu(x) G SO(N) for a.e. x G X. For instance, the map u(x) = x if X\ > 
and u{x) = (— x±, x 2 , ■ ■ ■ , Xjv) if xi < belongs to the Sobolev space W 1 ' 2 (Q,~R N ), 
Vu(x) G O(N) for a.e. x, but u is not an isometry. 

Since reformations have to preserve the volume, we have the following result. 

Theorem 3.8 (Rigidity). Let U C ~R N be an open connected bounded set. Let u : 
U — > M. N be a continuous, locally Lipschitz, open map such that C N {U) = C N (u(U)) 
and satisfying the following conditions 

(i) u(dU) C du{U) 

(ii) u is a.e. differentiable and Vw G 0(N) a.e. on U. 
Then u is an affine function. 

Proof. By (ii) it follows that u is locally a 1-Lipschitz function (see [14, Proposition 
3.4]). By the Area Formula and (ii) we infer 

C N (u(U)) = C N (U) = [ \Ju\dx= [ N(y,u,U)dy. 

Ju Jr n 

Therefore, it follows that N(y,u, U) = 1 for a.e. y G u{U). Observe that 

u(U) cR N \du(U) cR N \u(dU). (3.9) 

Therefore, for every x G U the topological degree deg(w(x), u, U) is well defined. Since 
u is a.e. differentiable, for a.e. x G U it results (see Lemma 5.9 of |19j ) 

|deg(tt(x), u, U) | = |sign (Ju(x))\ = 1. 

On the other hand, since u(U) is connected, by (13. 9p . u(U) is contained in a connected 
component of R N \ u(dU). Therefore, the degree is constant on u{U). It follows that 
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the sign of the Jacobian Ju is a.e. fixed. The conclusion follows by Liouville Theorem 
for Sobolev maps (see for instance [IS])- D 

For a related rigidity result involving local homeomorphisms see (59]. For quasi- 
isometries over Banach spaces see [91 Cor. 14.8]. 

Remark 3.9. Condition (i) holds of course for invertible maps u. If we are dealing 
with locally invertible maps, since continuous and locally invertible maps are open 
maps, actually by (i) the equality u(dU) = du{U) holds true. Moreover, if the map u : 
dU — > du(U) is injective, then u is globally invertible (see also [H]). Another classical 
condition for global invertibility holds for simply connected, or simply connectedly 
exhausted, target u{U) (see [59]). Moreover, suppose to have a continuous, locally 
invertible, surjective function u : X — » Y such that Vu(x) € O(N) for a.e. x. Then, 
D = N(y, u, U) is constant (see [3]) and by Area Formula we have 

DC N (Y)= [ N(y,u,U)dy= [ \Ju{x)\dx = C N {X). 

JY J X 

Hence, if C N (X) = £ N (Y), it follows that N(y,u,U) = 1 and hence u is globally 
invertible. 



By restricting to differentiable functions, if v is not isometric to /x then for any such 
reformation u it results ££{u) > 2. A natural question is if it is possible to get close 
to the minimum value by such reformations. 

We state in advance some stability properties for the inverse maps. 

Lemma 3.10. Let f n : X — > Y , f^ 1 : Y — )■ X be sequences of equi-Lipschitz mappings. 
U ||/n - /Hoc ->■ and H/- 1 - s-Hoo ->■ 0, then g = f' 1 . 

Proof. 

\g(f(x)) -x\< \g(f(x)) - f-\f(x))\ + - /- l (/„(x))| < 

< \\g-f- l \U + H\f{x)-f n {x)\ < Wg-f^ + HWf-U^ 
where if is a common Lipschitz constant for f~ l . Analogously we compute 

\f(g(y) -y\< \f(g(y)) - f n (g(y))\ + \f n (g(y)) - Wn\y))\ < 

< ||/ - / n |U + L\g(y) - f~\y)\ < \\f - / n |U + L\\g - f' 1 ^, 

for a common Lipschitz constant L for f n . Passing to the limit as n — > +oo we get 
the thesis. □ 

Lemma 3.11. Let f n : Q — > Y be a sequence of locally invertible mappings such 
that f n —>f locally uniformly on X. If the sequence (/~ 1 ) ng N is made of locally equi- 
Lipschitz mappings, then f is also locally invertible. 
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Proof. Fixed xo G X, suppose by contradiction to get two distinct sequences (x^hen, 
( x h)heN converging to x such that f(x\) = f{x\) \/h G N. Let e > be fixed. By 
uniform convergence, we find a large integer n such that \f n (x) — f{x)\ < e for every 
x G B(xo,r). As in jlTJ Proposition 7], we may assume that for every n G N f n is 
invertible in B(x ,r). Indeed, setting y = f(x ),y n = f n (x ), let U be the connected 
component of Xq contained in f~ 1 (B(yo,r)). By uniform convergence, we may assume 
that \ f n {x) — f(x)\ < | for every x G U . Since \y n — yo\ < |, there exists a connected 
component Q n of xo contained in /~ 1 (5(y ? y))- We claim that Q n C U. Indeed, if 
Xj —> x with Xj G Q n we have 

\f(x)-y \ < \f(x)-f n {x)\ + \f n {x)-y \ < ^ + |/„(x) - y \ < r. 

Therefore x G U. By local invertibility of f n , every path connecting y G B(y , ^) to 
may be lifted to a path in Q n starting at x . Hence f n (Q n ) = B(y , y). Since /„ 

is locally invertible in Q n , by Global Inversion Theorem ( [2] Theorem 1.8]), it follows 

that every f n maps homeomorphically Q n onto B(y , y). 

Let (7i be the connected component of x contained in / _1 (5(?/ , |)). For every 

X\ G U\ we have 

2r 

|/n(^i) - 2/o I < \f n (xi) - f{xi)\ + |/(ari) - 2/ | < y- 

Therefore, Z7i C /~ 1 (5(?/o, y)) which implies that Z7i C Q n . Therefore, the restriction 
of f n to U% is injective for every n > 1. Now, for a large ft, we may assume x\,x\ G 
5(x ,r) C Z7i. Then we compute 

\4 - 4\ = \fn\Uxl)) - fn\f n {xl))\ < H\f n {x\) - f n {xl)\ < 

H (|/n(4) " /(4)l + 1/(4) - /(4)l + 1/(4) - /n(4)l) < 2He, 
where if is a common Lipschitz constant for f~ l . By the arbitrariness of e we get the 
contradiction x\ = x\. □ 

Now we prove the following existence result for the variational problem of elastic 
reformation. 

Theorem 3.12. Let y, G V(9) and v G P(Y) so that y, = C N \_Q, v = C N L.Y . Then 
the variational problem 

minimize-f^w) | u G Ref(/i; v),u open s.t. HK < 2} (3.10) 
admits solutions whenever {u G Ref(/i; v),u open s.t. HK < 2} ^ 0. 

Proof. Since y = £ N \_Q, we may assume that X — O. Otherwise, we should assume 
that X fl B(x,r) is a convex, or a quasiconvex set. Let (w n ) ne N be a minimizing 
sequence. Given xq G Q, let K,H,r > as provided by Definition 12.21 It follows 
that the sequence (w n ) ne N is locally equi-Lipschitz on B(x ,r), see Lemma IL~2l or ^8] . 
Therefore, the sequence u n is pointwise equicontinuous on Q. By the Ascoli-Arzela 
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Theorem we extract a subsequence converging, uniformly on compact subsets of f2, 
to a continuous map u. For this continuous limit map u : Q — > ~R N it is easily 
seen that u#p = v. It remains to prove that actually u G Ref(/z;i/), namely that 
e u (x) < K(xq),c u (x) < H(xq) for every iGlfl B(xo,r). For X = Q, by Lemma [l~2l 
we get the Lipschitz condition 

\u n (xi) - u n (x 2 )\ < K(x )\x 1 - x 2 \ 

for every xi, x 2 G B(x , r) C VL. Passing to the limit asn-i +00 and then as r — > + , 
we obtain e u (x) < K(x ). Let X\ G B(x ,r) and y\ = u(x\) be fixed. Observe that 
by Theorem 12.171 the maps u n are locally invertible. Therefore, using Lemma [1.21 and 
Lemma [3.31 the inverse maps u~ l are also equi-Lipschitz. As in Lemma [3.111 we find 
the homeomorphisms u n : Q n — > B(yx,ri). By Lemma 13731 we get e u -i(y) < H(xq) 
for every y G B(yi,ri). Observe that u is open by the Domain Invariance Theorem. 
Hence u(Q) is actually an open set. By Lemma [1.21 it follows 

K 1 ^/) - u n l (Vi)\ < H(x )\y-yi\, G B(y 1 ,r 1 ). 

For a common neighborhood X\ G U\ C Q„ we have 

\x - x x \ < H(x )\u n (x) - M„(Xi)|. 

Passing to the limit as n — > +00 and then as x — > X\ we get c u [x-\) < H(xq). Therefore 
u G Ref (/i; v). 

Fixed e > 0, we find 5 > such that J E (H(x) + K(x)) dp < e whenever C N (E) < 5. 
By using a Vitali covering, we cover Q, up to a measurable set C N (E) = 5 > 0, by a 
finite number of disjoint neighborhoods £/j on which u n — >■ w uniformly and invertibility 
holds according to Lemma 13.101 and Lemma 13.111 
Since u#p = z/ we compute 



< Lip(u)(x)dp + / Lip(M _1 )(M(x)) d/i ) + / (H(x) + K(x)) dp 

i=l \JUi JU, J JE 

< 2_j\ / Lip(u)(x)dp + / Lip(-u _1 )(y) <iz/ J + e. 



i=i 

By Lemma [1.41 we get 



^(w) < / liminf ( / Lip(u n )(x)dp + / Lip(w n 1 )(?/) rfi/ I + e 

j-f n ^+°° XJUi Ju{U t ) J 

< liminf > I / Lip(it n )(x)<i/i + / Lip(u~ l )(y) dv I + e 
"^+°° ^ Vie/, A(^) / 

< liminf) / e Un (x) dp, + / c Un (x) dp )+£ = liminf 3&(u n ) 



+ e 
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Letting e — > + we get the thesis. □ 

Remark 3.13. In a similar way, according to Theorem 12.151 and Theorem l2.13l we obtain 
existence for the variational problem over the set {u G Ref (//; v), e n < \fi\ or the set 
{u G Ref(/x; v),u incompressible}. 

Remark 3.14. If the surjection property u(X) = Y is required for w G Ref(/x; z/), we 
may argue as follows. To check that u is onto, let us fix y G Y. Observe that u~ l 
are locally equi-Lipschitz. Arguing as in Lemma [3.111 for the sequence w" 1 we find a 
common neighborhood B(y ,r) such that 

u- 1 :B(y 0l r)^Uccn 

are simultaneously homeomorphisms. Therefore, since u n (X) = Y, we find a sequence 
x n — » x G U C Q such that M n (x„) = y . Then we have 

\u(x )-y Q \ < \u(x ) - u(x n )\ + \u(x n ) - u n (x n )\ + \u n (x n ) - y Q \ 

< \u(x ) - u(x n )\ + \\u- U n \\oo ->■ 

as n — > +00. 

Remark 3.15. Observe that the compactness of Ref(/i; z/) does not involve the energy 
functional M. In the setting of mappings with bounded distortion compare with [59l 
Ch. II Section 9]. In the case of H, K G L P (X) the above minimization result could 
be obtained by using Rellich-Kondrakov compactness in Sobolev spaces and the l.s.c 
of the p-Dirichlet energy. The point is that the above reasonings rely just on metric 
objects and can be generalized also in a metric framework by using the theory of 
Sobolev spaces over metric spaces (see Appendix IB"]) . 

Remark 3.16. For X compact, considering finite coverings, it turns out that H,K G 
L°°. Therefore, in such a case we may consider H, K as two universal constants. 
However the proof of Theorem 13 . 1 21 works as well for the non-compact case. It would be 
interesting to develop an analogous theory under weaker requirement on the functions 
H, K. 

We then characterize isometric measures by the following 

Theorem 3.17. Let n G V{Q) and v G V(Y), so that jj = £ N \_Q, v = £ N \_Y, for a 

given bounded set Y . Then, S(fi, u) — 2 if and only if there exists an isometry u such 
that u#fi = v. 

Proof. By Theorem 13 . 1 2 1 we get a minimizer u : Q — > R which belongs to Ref(/i; v). 
By Theorem 12 . 1 51 and Remark 12.151 it follows that u is globally invertible. By Lemma 
13.71 and Theorem 13.81 it follows that u is a local isometry, then (see for instance [9j Th. 
14.1]), u is an isometric map. 

□ 
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By Theorem 13.171 we have that by reforming a flat configuration \x in a corrugated 
one v it results S(fi, v) > 2. This last fact gives an alternative proof of the so-called 
Grinfeld instability (see [21]), indeed, by the changing of the geometry, any possible 
reformation must expand or contract some portion of the body. 

4. Generalized reformations 

The notion of reformation introduced in the previous section has some restrictions, 
indeed it is easy to show examples, like the one in Figure 4.1, in which every reformation 
has a large cost while allowing fractures of the body it is possible to map the initial 
measure by using local isometries. 

In 1 n 1 W 



ta h 



Figure 4.1. An isometric fractured reformation. 

Here we introduce a notion of generalized reformation. Our approach relies on 
measure theoretic tools mostly developed in the field of optimal mass transportation 
(see BQE21E3]). 

The notion of reformation map corresponds to the notion of the so-called transport 
map, i.e. u : X — > Y such that u#fi = v. A natural generalization of the transport 
map is given by the notion of transport plan. A transport plan between two probability 
measures \i G V{X) and v G V(Y) is a measure 7 G V{X x Y) such that 71^7 = 
71"^ 7 = u, where 7r*, i — 1, 2 denote the projections of X x Y on its factors. A transport 
map u corresponds to the the transport plan j u := (I x where I is the identity 

map of X. Observe that the set of transport plans with marginals \i and v, denoted 
by n(yU, u), is never empty since it always contains the transport plan 
We shall call generalized reformation, or reformation plan, of \x into v any transport 
plan 7 with marginals \x and v. 

Let us recall some known results which will play a crucial role in the following (we 
refer to 0H]). 

Definition 4.1. Let Ai(Y) be the space of Radon measures on Y . A map A : X — > 
M.{Y) is said to be Borel if for any open set B C Y the function x G X 1— > X X (B) is 
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a real valued Borel map. Equivalently, x i-> X x is a Borel map if for any Borel and 
bounded map if : X x Y — >■ R results that the map 



xEX^ J Lp(x,y)d\ x 



is Borel. 



Theorem 4.2 (Disintegration theorem). Let 7 G "P(X x Y) be given and let tt 1 : 
X x Y — > X 6e t/ie /irsi projection map of X x Y , we set fi = (tt 1 )^. Then for 
— a.e. x G X there exists v x G V(Y) such that 

(i) t/ie map x ^ v x is Borel, 

(ii) Vy> G C 6 (X x y) : J XxY ip(x,y)d-f = J x (J Y ip(x, y)du x {y)) dfi(x). 
Moreover the measures v x are uniquely determined up to a negligible set with respect 

to fi. 

Let 7 G Il(/i, v), as usual we will write 7 = v x ® //, assuming that v x satisfy the 
condition (i) and (ii) of Theorem 14.21 Obviously the transport plan \i <g) v corresponds 
to the constant map x 1— > v x = v. Let u : X — > Y, observe that for the transport plan 
7„ := (I x tt)#/x, the Disintegration Theorem yields 7^ = S u ^ <g) ji. 

Remark 4.3. Let X C K , we recall that the barycenter of a measure /i G V(X) is 
given by 

/9(/i) = x d/i. 



x 



If 7 = i> x ®fi, then, by Theorem 14.21 the map a; (->■ /3(z/z) is measurable. It is possible to 
define a generalized pointwise expansion and compression energy through the pointwise 
Lipschitz constant of the map x 1— >■ <£>(x) := /3(^ x )- Observe that for a transport map 
u, since /3(<y = x, we have 

r<p{x ) = r u (x ). 

However, it may happen that the map <p is an isometry although the target are quite 
far from being isometric as it happens in Figure 4.2 . 

In the sequel we will introduce the notion of generalized pointwise compression and 
expansion energy through the notion of 1-Wasserstein distance of measures. 

Definition 4.4. Let fi, v G V(X), the \-W asserstein distance between /1 and v is 
defined by 



W(ix,v)= inf / d(x,y) dj(x,y). (4.1) 

yeTl(fj,,u) J x 

Let us recall that by Kantorovich duality (see [U 1221 E2J E3]) the 1-Wasserstein 
distance between /i and v can be expressed as follows 



W{n, v) = sup <^ / tp d{n - u) I tp G Li Pl (X) } . (4.2) 



x 
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Figure 4.2. A barycenter isometric reformation. 



Lemma 4.5. The balls of (V(Y),W) are convex. 

Proof. Let fi G V(Y), r > be fixed, we consider z/ 1; v 2 G B := B(fi,r) C V(Y). For 
every i G [0, 1], let z/ 4 := tz/ x + (1 — t)z/ 2 . Then, by considering (14. 2p . for any fixed 
<p G Lip 1 (F), we compute 



ipd{y t -(j) = t ip d(ui - fi) + (1 - 1) \ (pd{v 2 — n) 
y Jy Jy 

< tW{v x ,Li) + {±-t)W{v 2 ,v)<r 

Passing to the supremum with respect to ip G Lip 1 (F) we get W(i>t, /■*) < r, hence 
v t eB\/te [0,1]. ... ^ 

The above lemma allows to pass from pointwise Lipschitz bounds to local Lipschitz 
condition over the metric space (V(Y), W) (see [15]). Let 7 = v x <g> fj,, the function 

f : X ^ (V(Y),W), f(x) = u x . (4.3) 

will be called disintegration map. Let us introduce the notion of generalized compres- 
sion and expansion energy in terms of the disintegration map /. 

Definition 4.6 (Generalized expansion and compression energy). For any reformation 
plan 7 = v x £g> fi of fi into v we define the pointwise expansion energy 

e 7 (x ) := limsup W ^^\ (4.4) 

x— ¥xq \X %0\ 

and the pointwise compression energy 

c 7 (x ) = limsup — ^ — ^-r. (4.5) 

By using H4.3[) we can state 

e 7 (x) = ej(x), c 7 (x) = Cf(x). (4.6) 
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The pointwise reformation energy is then defined by 

r 7 (x ) = e 7 (x ) + c 7 (x ). 

Remark 4.7. Notice that, since W(8 X , 5y) = \x — y\, if 7 is a reformation plan induced 
by a map u : X — > Y, say j u — (I x and / u is the disintegration map of 7, then 

it results 

r 7 (x ) = r /u (x ) = r u (x ). 

Definition 4.8. We define the set GRef(yU;z/) C n(/i, v) as the subset of reformation 
plans 7 of [i into v satisfying 

Vio G X : 3r>0,H,K s.t. e 7 (x) < K, c 7 (x) < H (4.7) 

/or every x G X fl -B(a;o, r). 

Remark 4.9. By ( I4.4p - fl4.6p the role played by the disintegration map is clear, hence 
one is led to argue as in the previous section trying to establish the analogous of 
Theorem 12.151 in the case of disintegration maps. Unfortunately in the general case 
of metric spaces some tools as degree theory are not available. Therefore, it is not 
clear if local invertibility follows by (14. 7\\ . Nevertheless, we shall see, analogously to 
Remark 12.151 that by restricting the analysis to the case of small reformations, as 
defined below, we can prove a global invertibility result suitable to the present case. 

Theorem 4.10. Let 7 G GRef(/x; v) and f : X — > V(Y) be the corespondent disinte- 
gration map. Assume HK < y/2. Then f is globally invertible on X. 

Proof. Following the proof of Theorem I2.15[ by the Metric Area Formula (see |6j H2J 
143] ) we have 

DU N {f{B)) = [ N(y,f,B) dU N {y) = [ J(MD(f,x)) dx < 

JV(Y) JB 

< [ e f (x) N dx < K N C N (B), 
Jb 

where for any seminorm P the metric Jacobian is defined by 
J(P) = Nu N ( [ P{v)- N dH N ~\v) 

Denote by g = f\ B . Since g is bi-Lipschitz and g(B) C f(B), we compute 

C N (B) = C N {g-\g{B)) < H N H N (g(B)) < H N U N (f(B)). 

Therefore, we get D < H N K N . Then the map / is globally invertible. □ 

Definition 4.11. Let us define the set of small reformation plans between \x and v as 
follows 

GRef (/i, v) = {7 G n(//, v) I e 7 < K, c 7 <H,HK< V2}. (4.8) 
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5. Finding reformation plans 

In the following examples we show that it is possible to compare shapes with regular 
disintegration maps despite no regular transport map does exist. 




B 



Figure 5.1. A disconnected target reformation 



Example 5.1. Consider a regular domain X C M. N splitted into Y = A U B for two 
disjoint regular domains A,Bc M. N in such a way 1 = C N (X) = C N (A) + C N (B). 
We find (see [311 [65]) two diffeomorphisms U\ : X — > A, u 2 : X — >■ B so that \ Ju\\ = 
C N (A), \Ju 2 \ = C N (B). 

Diffeomorphisms with constant Jacobian can be constructed by using the results of 
Indeed, let (p : Q — > Qi be a diffeomorphism. Assume for instance J(f(x) > 



Vi 6 fl and let f(x) = c N (Qi) J^p{ x )- Then 



f(x) dx 



% [ Jip(x) dx = C N (Q). 
i) Jq 



By the results of [H], there exists a diffeomorphism u : — > such that Jw 
Setting ip — (p o u^ 1 : O — >- it follows that J-0 = 75 n mk ■ 

Let = C N (A)S Ul ( x ) + C N \B)8 U2 (x)i then the reformation plan 7 := v x <g> ji has 
(i = £ N \_X and z/ = C N L.Y as marginals. We claim that the function f(x) = v x is, 
at least locally, bi-Lipschitz. Indeed it results 

W(y x , u X0 ) = C N (A)\u\{x) - Ul (x )\ + C N (B)\u 2 (x) - u 2 (x )\. 
Since u\, u 2 are diffeomorphisms, we find constants K\ 2 , H\ y2) H, K such that 



1 . 



X \ 



< 



C N (A), 



\x 



C N (B) 
Xq\ H — — \x 



x \ 



H\ H 2 
< C N (A)\ Ul (x) - Ul (x )\ + C N (B)\u 2 (x) - u 2 (x )\ 
= W(v x ,v Xo ), 



28 



LUCA GRANIERI, FRANCESCO MADDALENA 



W(v x ,v X0 ) = C N (A)\u 1 (x)-u 1 (x )\ + \C N (B)\u 2 (x)-u 2 (xo)\ 

< C N (A)Kx\x — xq\ + C N (B)K 2 \x — x \ 

< K\x — xo\. 

Remark 5.2. The above construction is possible also for a class of star-shaped domains 
as in [20l Theorem 5.4] by considering bi-Lipschitz maps in place of diffeomorphisms. 

Moreover, generalized reformation maps are useful to compare near-isometric shapes. 





1 

n 







Figure 5.2. Bending a rectangle. 

Example 5.3. Consider a rectangle R and a bended one with the bended size of K 
Consider the maps 

uUx) = ( 1 - - ] (Ax + a) , u 2 (x) = - (Bx + b) 
\ n J n 

for orthogonal matrices A, B and then the reformation plan 

7 = (( X ~ n) 6ui(x] + ^ U2{3;) ) ® /i ' 
where /x = L i?. We compute 

W(u x , v xo ) = ( 1 - - J |Mi(x)-Mi(x )|+-|M 2 (a;)-M2(a;o)| = ( ( (1 - - J + -4 | |x-x 
\ n J n \\n/n z l 

Therefore the function /(#) = ^ x is, at least locally, bi-Lipschitz and 

e 7 (x ) = (1 - - J + ^ -)> 1 
\ n J n z 



as n — > +00, while c^(xq) 
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e 7 (xo) ' 



H 1 h 



y = x + 1 



2^2 



i 1 h 



y = —x — 1 



Figure 5.3. An horizontal segment, with mass 2 a/2, splitted into two 
different ones. 



Example 5.4. Consider the situation displayed in Figure 5.3. Defining v x — \ (<Wi + 8-x-i) 
we find W(v x , u Xo ) = \/2\x — xo\. Hence r 7 = e 7 + c 7 = y/2 + 

Example 5.5. Let X C be a measurable set with C N (dX) = 0. We find an 
increasing sequence of polyhedral sets X n such that X = U n>1 X n up to a negligible 
set. Let Y C M. N be the unitary cube, C N (X) = C N {Y) and let Y n C F be a rectangle 
such that C N (Y n ) = C N (X n ) Vn G N. Let /1 = C N LX, v = £ N \_Y. 
We find a sequence ( - u n )neN so that Vn G N «„ : X n — )■ F n is a bi-Lipschitz map with 
Ju n = 1. The volume constraint implies that X n := Lip(w n ) < fT, if n := Lip(w^ 1 ) < 
H. In particular, for every x,y G X„ we have 



1 



k - V\ < l«n(a?) - u n (y)\ < K\x - y| 



By Lipschitz extension, we may consider u n as defined on the whole X. By Ascoli- 
Arzela Theorem we find m„ m uniformly. It follows that 



H 



\x — y\ < — u{y)\ < K\x — y\, 
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up to a zero measure set. Moreover 



f(u(x)) dx = lim / f(u n {x)) dx 



x lt — ? ~l~ OCJ j x 



lim ( / f(u n (x))dx+ / f(u n (x))dx 

X n J X\X n 



lim / f(y) dy= f(y) dy. 

'Y„ JY 



Hence, = v. 

6. Variational problems on generalized reformations 

The notion of generalized reformation involves the Lipschitz pointwise constant of 
maps in a metric space framework. For the associated integral energies it is natural to 
consider some notion of Sobolev spaces in a metric setting. There exist different notions 
of such metric Sobolev spaces which coincide provided some mild assumptions such as a 
doubling condition, a Poincare inequality and a power of integrability 1 < p < +oo are 
satisfied. We refer the reader to the Appendix B and the references therein for further 
informations. In particular the requirement on the power 1 < p will be important to 
state a general existence result, see Theorem 16. 8 \ for the variational problem related to 
generalized reformations. Actually, these kind of assumptions seem to form a natural 
context to work with in the setting of metric analysis. Therefore, along all this section 
we will assume 

X = tt cR N compact and satisfying flR3l) and flBTil) . 



(6-1) 

Y C M. N compact. 

Definition 6.1. Let 7 G n(yU, v). We define the reformation energy of ^ as follows 

^( 7 ) = / (c 7 + e 7 ) dfi. (6.2) 



J x 

Remark 6.2. With a slight abuse of notation we are using the same symbol & to 
denote the reformation energy functional defined on the space of reformation maps 
and the analogous defined on the space of reformation plans. Since in the paper it 
always appear with its argument specified, there is no risk of confusion. 

Theorem 6.3. Let 7 = f(x) ®/i G GRef(/z; v) be such that = 2, \i absolutely 

continuous with respect to the Lebesgue measure. Then there exists an open dense 
subset of X on which the disintegration map f is a local isometry (with respect to the 
Wasserstein distance). 

Proof. First observe that since X is quasiconvex (see for instance jl5l Lemma 6.1]), 
then / is a Lipschitz function. We have e 7 = c 7 = 1 a.e. By J5TJ Prop. 1.1, Sec. 3], 
there exists an open dense subset U C X on which / is locally bi-Lipschitz. Therefore, 
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consider a bi-Lipschitz map / : B — > V(Y) for an open ball B C U. For x%, x 2 G B, by 
using Fubini Theorem, we find a curve rj connecting Xi,x 2 as in [TU Prop. 3.4] in such 
a way for a.e. t it results e 7 (^(t)) = 1 and £(77) < |xi — x 2 | +£. Since / is Lipschitz, the 
curve p : [0, 1] — > (V(Y), W), defined by p t = f(rj(t)) is Lipschitz too. Hence, it admits 
a tangent vector v (see Theorem IA.2I) . Fixed u G Lip 1 (F), by standard approximation 
argument we may suppose that u G C 1 . Therefore, by using the continuity equation 
(1A.4I) we compute 

u d {fixt) - f{x 2 )) = J u d{p x - po) = J ^(J u dpt^j dt = 
/ (du,v)dp t dt< / / \v\dp t dt = / \p\(t) dt < / ey(r} t )\f]\dt < 

>Y JO JY JO JO 

< l(rj) < \xi - x 2 \ + e. 
Taking the supremum with respect to u and letting e — > + we get the 1-Lipschitz 
property 

W(f(x 1 )J(x 2 ))<\x 1 -x 2 \. 

To get the opposite inequality, we argue as follows. Set po = f(xi),pi = f(x 2 ), let us 
consider a geodesic p t : [0, 1] — > V(Y) between p and p±, i.e. l(p) = W(f(xi), f(x2))- 
Since / is bi-Lipschitz, there exists an injective Lipschitz curve 7 : [0, 1] — > B con- 
necting Xi,x 2 such that p t = f(j(t)). Again by using a Fubini type argument, we 
find a sequence of Lipschitz injective curves (7„)„ e N so that j n —> 7 uniformly and 
Lip(/~ 1 )(/(7n(^))) = 1 for a.e. t G [0, 1]. Therefore, we get a n = /(7 n ) — > p uniformly 
in (V(Y),W). By the upper semicontinuity of the Hausdorff measure along the se- 
quence a n (see for instance [101 Lemma 4.1]), recalling that for injective curves it results 
1(a) = H 1 (cr([0, 1])) (see [7]), fixed e > 0, we find a Lipschitz curve a connecting p Q 
and pi such that Lip(/~ 1 )( ( r(t)) = 1 for a.e. t G [0,1] and /(a) < W{f{x 1 ),f(x 2 ))+e. 
Finally, we compute 

\x 1 -x 2 \ = \f-\a(0))-f-\a(l))\= r <! 



dt f-\°{t)) dt < 
< [ \&\ w (t)dt = l(a)<W(f(x 1 ),f(x 2 ))+£. 

JO 

Letting e — > + we get the thesis. □ 

Theorem 16.31 should be compared with Theorem 3.17. The main restriction is on 
invertibility which is just on an open dense subset. We may say that this open set 
is of full measure, actually coinciding with the whole space, just for the case of small 
reformations as done in Theorem 16.41 below. There are different restrictions in doing 
so for the general case. A first matter relies in characterizing the set where a map is 
locally invertible on a metric setting. A second one relies on the fact that the integral 
functional M gives a.e. informations, while invertibility requires global conditions. 
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Therefore the matter is on passing from a.e. conditions to everywhere ones. In the 
results concerning transport maps, this difficulty was overcome by using degree theory 
in Wt N . Therefore, something similar to degree theory over metric spaces should be 
needed in order to handle with this kind of questions. 
Let us introduce the notation 

£ G (ji, v) = inf{^( 7 ) | 7 G GRef (/z; u)}. (6.3) 

Concerning symmetry properties of the above generalized reformation energy, the same 
reasonings made for transport maps, compare with cerning Definition 13. 4[ hold as well. 
We remark here that this time the question of symmetry is not just a question on 
invertibility. For instance, the transport plan 7 = / <g> \x between \x and u, considered 
in figure 4.1 is isometric, i.e. W(f(x), f(x )) = \x — x \. However, reversing the target 
measures we see that the transport plan between v and fi is just locally isometric and 
no transport plan g®v between v and \x is isometric. The fact is that the corresponding 
disintegration maps are of the form 



g:Y^V(X). 

Therefore, symmetry questions are quite involved and here we do not further consider 
them. 

We state the following characterization of the lowest possible value of the generalized 
reformation energy. 

Theorem 6.4. If Sg{^, v) = 2, with fi absolutely continuous with respect to the 
Lebesgue measure, then the infimum is attained at a local isometric reformation plan. 

Proof. Let 7 n be a minimizing sequence. By compactness of V[X x Y), by passing to 
a subsequence, we may assume that 7„ — ^ 7. It follows that 7 is also a transport plan 
between \i and v. By disintegration, we also assume that 7n = f n (x) <8> fx, 7 = u x <g> \i. 
For any fixed (p G C(X), ip G C(Y), we get 



<p(x) / ip(y)dv x dn= (p(x)ip(y) d>y = lim / ip(x)ip(y) d^ n = 

= lim / tp(x) ( [ i>(y)df n (x)) dfx. (6.4) 

n^+OO J X \Jy J 

By density of continuous functions, it follows that f Y ip(y)df n (x) — ^ f Y ip(y)di/ x in 
Lebesgue spaces of integrable functions. 

Since X is quasiconvex, by definition of generalized reformations, it follows that 
the sequence f n is equi-Lipschitz on X. By Ascoli-Arzela Theorem, by passing to 
a subsequence we have that f n —>f uniformly. Since the disintegration is uniquely 
determined, it follows that f(x) = v x for Indeed, since the Wasserstein 

distance metrizes the weak* convergence of measures (Y is compact), for every ip G 
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C(Y) we have 

ipdf n {x) -> / 4>df{x). 



Hence, for every (p G C(X), passing to the limit under the integral sign and by (16. 4p 

we get 



tp(x) yj ipdf(x)J d^ = ^lim^y tp(x) ipdf n (x) J dp = J <p{x) yj ipdv x J dfi. 

Observe that each /„ is locally invertible on an open dense U n C X. Since we are 
dealing with small reformations, by Theorem 14.101 it follows that each /„ is globally 
invertible and locally bi-Lipschitz on X by Lemma l4~5l By Lemma \3.11\ we have that 
also / is locally bi-Lipschitz on X. Moreover, it easily seen that / G GRefo(/U;^). 
Since 



2< / (e 7n + — U/i<^( 7n ) VnGN, 
passing to the limit we get 

lim / g n {x)dfi = 2, 

n— >+oo J 

where g n (x) = e 7n + Passing to a subsequence we have g n — > 2 a.e. Since 

9n( x ) — f&n) f° r Pit) — t + t, by continuity of cp it follows that e 7n — > 1 a.e. On 
the other hand, c 7 > — yielding lim inf n _s. +00 c 7 > 1 a.e. Since -y n is a minimizing 
sequence for we get 

2 = lim &("y n ) = 1 + lim / c ln dfx. 

n— >+oo n— >+oo J ^ 

and by Fatou Lemma we infer 



1 < / lim inf c 7n dfi < lim / c 7n = 1. 

Therefore, by passing to a subsequence, we also have that c 7n — > 1 a.e. Arguing as 
in the proof of Theorem 16.31 we locally find in Q a curve rj : [0, 1] — > V(Y) such that 
e ln (r](t)) — > 1 a.e. and l(rj) < \x\ — x 2 \ + £■ Therefore we get 

W(/n(zi),/n(* 2 )) < [ ^n(v(tm\(t)dt. 
JO 

Passing to the limit we obtain 

W(f(x 1 )J(x 2 ))<l( V ) < \ Xl -x 2 \+e. 
Letting e — > + we obtain the 1-Lipschitz condition 

W(f{x 1 ),f{x 2 )) < \ Xl -x 2 \. 
Arguing again as in the proof of Theorem 16.31 we locally obtain 

W{f{xi)J{x 2 )) = \X!-X 2 \, 
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hence = 2. 



□ 



Remark 6.5. To recover a global isometry in the above results as in Theorem 13.171 
one should establishes some metric version of Liouville Rigidity Theorems as in The- 



A natural question concerns the validity of an existence result as in Theorem 13.121 
However, we observe that the approach pursued in the proof of such result involve the 
push-forward of the transport map. Therefore, for generalized reformations, the push- 
forward of the disintegrations maps is involved. This point of view leads to consider 
a variational problem over transport classes as introduced in [3D]. The definition of 
transport classes is the following 

Definition 6.6. Let 7, rj G n(/x, v) with 7 = f(x) ®fi,r] = g(x) <8> \i be given. We shall 
say that 7 and 77 are equivalent (by disintegration), in symbols 7 ~ 77, if /#// = g#fi. 
For any 77 G n(/j, v) with 77 = g(x) <8> fi, we shall call transport class any equivalence 
class of a transport plan 77 and it will be denoted by [77], i.e. 



For a transport map u the disintegration map is given by 1 ^ S u m. In [SUj it is 
shown that every such disintegration map leads to the same push-forwarded measure. 
In other words, all the reformation plans of the form (I x u)#/J, belong to the same 
transport class. Moreover, the following result holds true 

Proposition 6.7. Let u : X — >■ Y be such that u#fi = v and let 77 = (I x = 

8u(x) <H> t 1 - If 7 € [v] then there exists v : X — > Y such that 7 = 5 V ( X ) <8> fi, i.e. the 
transport plan 7 is concentrated on the graph of v. 

In this perspective, the variational problem A3. 10[) studied in section 3 could rewrit- 
ten as 

minimize{^ , (u) | u G Ref(/x; v)} = minimize G Ref( At ; ! /){^ , (7) I 7 ^ [8u(x) ® lA}- (6.6) 

However, by passing to transport plans, different transport classes arise. By the above 
discussion it seems natural to fix a transport plan rj G n(/i, u), r\ — g{x) (g) \i and to 
consider the variational problem 



Theorem 6.8. (Existence of optimal reformation plans) Assume (16. ip and let 77 G 
GRefo(/x; v) be given. Then, for every p > 1 the variational problem 



orem l3T8l 



W = {7 = f(x) ® /i I /#// = 



(6.5) 



minimizeGR e f(^;^) {^(7) I 7 e [77]} . 



(6.7) 




(6.8) 



admits solutions. 
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Proof. Let 7„ = f n {x) ®/i be a minimizing sequence. Let f n — >■ / uniformly with 
respect to the Wasserstein distance as in the proof of Theorem 16.41 By Lemma 11.41 we 
get the lower semicontinuity of the term J x e^(x)dpi. Moreover, by Lemma [3.31 we get 

/ <*{x)d l i= [ Lip p (r 1 )(/(x)) rf/i (6.9) 
Jx Jx 

Since (16. ip X satisfies the doubling condition given in Definition IB. 41 and the Poincare 
inequality given in Definition IB.5[ we can apply the theory of Sobolev spaces over 
the subset f(X) of the metric space (V(Y), W, f#p) (see Appendix B). Moreover (see 
|56j). since for p > 1 the pointwise Lipschitz constant Lip(g) is the minimal generalized 
upper gradient of the locally Lipschitz map g ([551 Theorem 5.9]) and the Cheeger p- 
energy (IB. ip is lower semicontinuous with respect to LP convergence ([23 Theorem 
2.8]), by using (I6.9P we have 

f t*(x)dn = [ U^if-^y) d(M < liminf [ Lip^/" 1 )^) d(/ #/ i). (6.10) 
By taking into account the condition (f n )#p = f#fi Wn G N, we get 




□ 



Appendix A. Curves in metric spaces 

For reader convenience here we just summarize some basic results. For analysis in 
metric spaces we refer to [71 [3H [35] . For Lipschitz function on a metric space (M, d) 
we introduce the metric derivative according to the following definition. 

Definition A.l. Given a curve p : [a, b] — > (X,d) we define the metric derivative at 
the point t <E]a,b[ as the limit 

h^O h 

whenever it exists and in this case we denote it by \p\{t). 

Of course, the above notion of metric derivative coincides with the metric differential 
(11.21) . If p : [a, b] — > (X, d) is a Lipschitz curve, by metric Rademacher Theorem the 
metric derivative of p exists at £ 1 -a.e. point in [a, b]. Furthermore, the length of the 
Lipschitz curve p is given by 

l{p)= f\p\{t)dt. (A.2) 



We restrict to the case of V(X) := (V(Q),W). The following theorem relates abso- 
lutely continuous curves in V(X) to the continuity equation. 
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Theorem A. 2. Lett H- pt G V(X),t G [0, 1], be a curve. If p t is absolutely continuous 
and \p\ G L 1 (0, 1) is its metric derivative, then there exists a Borel vector field v : 
(t,x) i — y v t (x) such that 

v t G LP{X,p t ) and \\v t \\ LP{x , pt) < \p\(t) for C 1 - a.e. t G [0,1] (A.3) 

and the continuity equation 

p t + div(vp t ) = in (0, 1) x X, (A.4) 

where the divergence operator is understood with respect to the spatial variables, is 
satisfied in the sense of distributions. 

Conversely, if p t satisfies the continuity equation (1A.4j) for some vector fields v t such 
that \\Vt\\LP{p t ) G -^ 1 (0, 1), then t i— >■ p t is absolutely continuous and 

\p\{t) < \\v t \\Lp(x,p t ) for C 1 - a.e. t G [0, 1]. 

Remark A.3. The minimality property fl A. 3j) uniquely determines a tangent field v t . 
We will refer to v t as the tangent vector associated to the curve t >-)■ p t . The continuity 
equation flA.4j) has been used in the Monge-Kantorovich theory since its beginning for 
many applications. The fact that it characterizes the absolutely continuous curves 
on the space of probability measures equipped with the Wasserstein metric was only 
recently pointed out and the full proof is contained in [5]. 

An immediate consequence of the continuity equation is the following 



Lemma A.4. For every solution (pt,Vt) of the continuity equation ( A.4 ) and for every 
f G C X (X) it results 

d 



(jf yj f(x)dp t J=j (d x f(x),v t (x))dp t (A.5) 

in the sense of distributions. 

Actually, it turns out that the map / m- j x fdp t belongs to W{£(0, 1). Therefore, 
formula (1A.5j) holds for a.e. t G (0, 1). We refer the reader to [U [29] for proofs and 
more details. 

Appendix B. Sobolev spaces on metric spaces 

There are several ways to generalize the notion of Sobolev spaces into a metric 
framework, see for instance [121 HI H21 |36j H5j |56j |6T]. The approach based on the 
notion of upper gradient (see [121 ESI ESJ |6T]) seems to be more appropriate to the 
context of this paper. 

Definition B.l. Let (X,dx), (Y,dy) be metric spaces, let U C X be an open subset 
and let u : U —*Y be a given map. A Borel function g :[/—>- [0, +00] is said to be an 
upper gradient of u if for any unit speed curve 7 : [0, /] — > X it results 

dy(«(7(0)),u(7(0))< I 9{l{s)) ds. 
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If u : U — > Y is Lipschitz, then the pointwise Lipschitz constant Lip(w) is an upper 
gradient for u, see [121 [321 EI]- For w G £ p (t/, V"), the Cheeger type p-energy is defined 
as follows 

E p (u) = inf liminf |0 n |£ p , (B.l) 

where the infimum is taken over the sequences (u n , g n ) such that g n is an upper gradient 
of u n and u n — >■ w, g n — >■ (7 in L p . By definition ( IB. II) it immediately follows 

£L(w) < lim inf E p (u n ) (B.2) 

whenever u n — > u in L p . The Cheeger metric (l,p)-Sobolev space is defined as 

H lj, (U, Y) = {u G L P {U, Y) : E p {u) < +00}. 
We need two more definitions. 

Definition B.2. A function g G L p is called a generalized upper gradient for u G 
H 1,P (U,Y) if there exists a sequence (u n ,g n ) such that g n is an upper gradient for u n 
and u n —> u, g n —> g in LP . 

From Definition IB. II it follows that \g^ LV > E p (u) whenever g is a generalized upper 
gradient for u. 

Definition B.3. A generalized upper gradient g for a map u G H 1,P (U,Y) is said to 
be minimal if it satisfies \g\ p LP = E p (u) 

Under some regularity requirement on the target metric space Y, it may be proved 
(see [56]) that every u G H l > p (U,Y), with 1 < p < +00 admits a unique minimal 
generalized upper gradient g u . This minimal generalized upper gradient coincides 
with the pointwise Lipschitz constant Lip(w) under some geometrical property of the 
spaces (X,n),Y (see [561 Theorem 5.9]). In particular, a crucial role is played by the 
doubling condition and a weak Poincare (l,p) -inequality for the space (X,/i). 

Definition B.4. A measure [i over a metric space X is said to be "doubling" if fi is 
finite on bounded sets and there exists a constant C such that for every x G X and 
every r > the following inequality holds 

v(B(x, 2r)) < Cfx(B(x, r)). (B.3) 

Definition B.5. Let 1 < p < +00. A metric measure space (X,d,fi) is said to satisfy 
the weak Poincare (l,p) -inequality if, for any s > 0, there exist constants C, A > 1 
such that, for any open ball B(x,r) with < r < s, function f G L 1 (5(x, Ar))) and 
upper gradient g : B(x,Ar)) — > [0, +00] for f , the following inequality holds 

dfi<C[-f g p d/i) (B.4) 



B(x,r) 



f 



f dfi 



B(x,r) 
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Observe that under some geometrical requirement on X, the Poincare inequality 
(IB.4I) may be required to hold just for Lipschitz functions / (see [551 ES])- The eu- 



clidean space M. N equipped with the Lebesgue measure C N is doubling and satisfies 
the above Poincare inequality with A = 1. Given a square Q and fi = £ N \_Q, by the 
inequality 

^£ N (B(x,r)) < n{B{x,r)) < C N (B(x,r)), 

holding for every ball B(x, r) of Q and the usual Poincare inequality on convex sets, it 
follows that (Q,n) is doubling and supports the Poincare inequality (1B.4[) . Since the 
doubling condition and the Poincare inequality are stable under bi-Lipschitz maps, 
every diffeomorphic (or bi-Lipschitz), with volume preserving maps, domain Q (as 
balls, see for instance [20J [31]) with the same volume of the square Q, equipped with 
the measure v = C N Lfi is doubling and supports the Poincare inequality flB.4j) . For 
more details on the doubling and Poincare inequality we refer the reader for instance 

to mug una as]. 
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